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Abstract 

We show that Af — (2, 2) SU(N) super Yang-Mills theory on lattice does not have 
sign problem in the continuum limit, that is, under the phase-quenched simulation 
phase of the determinant localizes to 1 and hence the phase-quench approximation 
becomes exact. Among several formulations, we study models by Cohen-Kaplan- 
Katz-Unsal (CKKU) and by Sugino. We confirm that the sign problem is absent 
in both models and that they converge to the identical continuum limit without fine 
tuning. We provide a simple explanation why previous works by other authors, which 
claim an existence of the sign problem, do not capture the continuum physics. 
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1 Introduction 



Super symmetric Yang-Mills theory (SYM) attracts broad interests as a candidate of the 
physics beyond the standard model pQ. It is also a promising candidate for the nonper- 
turbative formulation of the superstring theory [21 HI El E] • Although it is important to 
study its nonperturbative properties for both applications, and one of the most successful 
nonperturbative approaches for gauge theory is the lattice simulation, however, notorious 
difficulty of the lattice SYM prevented it for long time. Sometime ago there appeared a 
breakthrough for two-dimensional theories of extended supersymmetries [H El [TUl ITT] , 
(see also [12] for a recent review); the correct continuum limit is obtained without param- 
eter fine tunings for most of these models, at least at perturbative levelQ Furthermore 
it turned out that by combining two-dimensional lattice with matrix model technique 
(fuzzy sphere) |14j , 4d J\f = 4 theory is constructed without relying on the parameter fine 
tunings to all order in perturbation theory |15j . (For other very elegant formulation in 
the planar limit, which preserves 16 supersymmetries manifestly, see [16]. This method 
is applicable to various kind of theories for which lattice and/or fuzzy sphere technique 
are not applicable [E].) These formulations provide robust ways to test the AdS^/CFT^ 
correspondence, the gauge/gravity duality JSJ E] and matrix string conjecture |4] at non- 
perturbative level. In one dimension, such a program has been pursued extensively by 
using a non- lattice method [18] and the gauge/gravity duality E] has been confirmed 
very precisely [191 E01 EU [221 E3] , including the stringy a' corrections [211 E3] • The lattice 
simulations are also applicable for this system and qualitatively consistent results have 
been obtained [24] . We can expect study of two and four dimensional theories will pro- 
vide even richer insights; in two dimensions simulations in this context are already ongoing 
|25 l l26 l l2"7] and it is urgent to establish the validity of the lattice models at nonperturbative 
level by detailed simulations. 

However there is a possible obstacle for simulations: supersymmetric theories with 
eight and sixteen supersymmetries suffer from the sign problem [28J1H. On the other hand, 
with four supersymmetries (i.e. 4d N = 1 pure SYM and its dimensional reductions), there 
is no sign problem. It can easily be seen as follows. In Weyl notation, with an appropriate 
choice of the gamma matrices, the Dirac operator is written as 



where a = —H2 and a l (i = 1,2,3) are Pauli matrices. By using a 2 {ia^)a 2 = (ia^ 1 )* and 
the fact that is real in adjoint representation, we obtain 



1 There is a lattice formulation of 3d maximally supersymmetric Yang-Mills without fine tuning [10] , 
For 4d TV = 1 pure supersymmetric Yang-Mills, the chiral symmetry guarantees supersymmetric continuum 
limit. Recent results are found in [13] . 

2 In the case of maximally supersymmetric matrix quantum mechanics, agreement with the dual gravity 
prescription has been observed by ignoring the phase of the Pfafnan, even when the sign fluctuates violently 
20, 22]. It has also been observed that the PfafEan is almost real positive for 5(7(2) [23] . It would be nice 

to understand why it happens, but it is out of scope of the present paper. 

3 Even if there is a sign problem, measurement of the sign factor itself is interesting because it is related 
to the Witten index with a suitable normalization |29| . 



D = ia^D 




a 2 fta 2 = D*. 



(2) 
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Therefore, if ip is an eigenvector corresponding to an eigenvalue A, a 2 ip* is also an eigen- 
vector, with eigenvalue A*. They are linearly independent and eigenvalues appear in a pair 
(A, A*). This assures the positivity of the determinant after removing A = modes. 

At discretized level, positivity of the determinant can be lost. In zero dimension (ma- 
trix model) |30j . there is no sign problem, because there is no need for the regularization. 
In Id theory and fuzzy sphere construction of three and four dimensional theories [16) . by 
using the momentum cutoff prescription [18] the determinant is positive even at discretized 
level [31]. 

In lattice constructions of two-dimensional SYM, the determinant is in general complex 
at discretized level. In Sugino's model [8] for two-dimensional theory, however, the sign 
problem disappears as one approaches to the continuum [321 S3]; that is, if one performs 
the phase-quenched simulation, distribution of the phase factor of the determinant)^ peaks 
at 1 in the continuum limit. Therefore the phase quench approximation becomes exact at 
continuum. In this case, in addition to the absence of sign, the agreements with analytic 
calculations in small volume region have been observed as well [25J by using techniques 
developed in [34J . Therefore the absence of the sign is the property of the correct continuum 
limit, as expected. Numerical studies of this model can also be found in [35) 136]. 

On the other hand, for other two-dimensional lattice models which are supposed to 
have the same continuum limit, an existence of the sign problem has been reported |37) 
138] . In Cohen-Kaplan-Katz-Unsal (CKKU) model, which is equivalent to a model by 
D'Adda et.al. [llj with a specific choice of parameters, Giedt reported that the determinant 
has a complex phase, and the phase fluctuates violently if one chooses random lattice 
configurations [37] , However the importance sampling has not been performed in [37J, 
and hence this result has nothing to do with the continuum limit as the author remarked 
correctly. This model was studied later with importance sampling in [38^, where the 
sign problem was reported as well. However, it is not clear whether it is a property of 
the continuum, because they could not evaluate physical quantities because of the "sign 
problem" and hence could not estimate how close to the continuum limit they have reached. 

In this paper, we resolve the confusion mentioned above. We show the absence of sign 
problem in the Sugino model and the CKKU model in the continuum limit. (In Sugino 
model, the absence of sign has been explicitly reported for SU{2) theory in [32) 133] . For 
N > 2 theory, we briefly checked but have not mentioned it in [25], because our emphasis 
was put on other physical quantities. In this paper we show the detail for SU(N) with 
N > 2, together with new data for SU(2).) The action in the continuum is obtained from 
4d M = 1 SYM through the dimensional reduction, and is given by 

S = - J dxj^ dy Trl-F%, + -{D^f - -[X^Xtf - -$T»D^ - -^[X.M , 

(3) 

where fi and v run x and y, i and j run 1 and 2, and T 1 = (r^,r*) are gamma matrices 
in four dimensions. Xj are N x N hermitian matrices, ip a are N x N fermionic matrices 

4 Strictly speaking, in Sugino model fermions are Majorana and hence we calculate the Pfafnan. 

5 In [38], the bosonic fields are denned as e z where z is a complex field, and thus a natural way of 
extracting physical quantities is different from the original CKKU model. For SU(N) gauge group as in 
|38| . this definition is different from a general complex field originally defined in [TJ. We would like to 
thank S. Catterall for detailed explanations on his work. 
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with a Majorana index a and the covariant derivative is given by = —i[Au, ■ ]■ The 
only parameters of the model are the size of circles L x and L y . (Note that the coupling 
constant can be absorbed by redefining the fields and coordinates. Therefore we take the 
't Hooft coupling A to be 1. Then the strong coupling corresponds to the large volume.) 
We study this system by using two lattice models (CKKU and Sugino) numerically and 
show the absence of the sign problem. We evaluate expectation values of some physical 
quantities and see that the results show nice agreements. Note that small volume behavior 
of the Sugino model is consistent with known analytic estimates [25J . 

One obstacle for the simulation is the existence of the flat direction, along which 
two scalar fields X\ and X2 commute. In contrary to a theory on R 1 ' 3 , there is no 
superselection of the moduli parameter in this case. That is, eigenvalues of scalars are 
determined dynamically. Therefore, some mechanism which restrict eigenvalues to a finite 
distribution is necessary for the stable simulation. In addition, to obtain an interesting 
dynamical system, having a (small) finite region for the eigenvalues is important as well; if 
the eigenvalues of the scalar spread so large, the theory would run into the abelian phase, 
which is just a free theorj^. In this work, we introduce soft SUSY-breaking mass to scalar 
fields 



so that the flat direction is lifted_|. It is crucial to control the flat direction for various 
reasons. We have just mentioned two of them — stability of the simulation and interest- 
ing non-abelian phase. There is one more; in order to guarantee the correct continuum 
limit, the eigenvalue must be smaller than the cut off scale ~ 1/a. Especially in the 
CKKU model, we can decompose the bosonic field to appear scalars X as a log of positive 
Hermitian variables H |39|, 



where o is the lattice spacing. In order to obtain the tree level action, one has to assume 
aX <C 1, expand it in powers of aX and neglect higher order terms. Therefore, unless 
aX <C 1, one cannot get to the continuum limit even at tree level. Actually it turns out 
that one of previous works, whose conclusion contradicts with ours, does not satisfy this 
condition. 

This paper is organized as follows. In § [2] we study the CKKU model. We introduce 
the model in § 12.11 and then consider a structure of the Dirac operator in § 12.21 Then we 
show the absence of the sign problem in § 12.31 In § [3] we show the absence of the sign 
problem in Sugino model. Then in §2] we confirm that two models (CKKU and Sugino) 
converge to the same continuum limit. To our best knowledge, this is the first result from 
the CKKU model in the continuum limit and shows in fact we can take a supersymmetric 
continuum limit without any fine tunings. In §[5] we explain why previous works by other 
authors fail to capture the continuum physics. 

6 Which phase is preferred is in fact a dynamical question. At large- N, the fiat direction is lifted and 
the system stays non-abelian phase; see [25] . This phase is an analogue of the black 1-brane solution in 
type IIB supergravity. 

7 For 8 and 16 SUSY models, there exists SUSY-preserving mass deformation [15] . 




i=l,2 



(4) 




H = exp(aX) 



(5) 
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2 CKKU model 



In this section we study the CKKU model. In § !2.1l we introduce the model. Then in § 12.21 
we discuss the structure of light modes in the model, which is crucial for the analysis of 
the sign problem shown in § 12.31 



2.1 The model 



Here we consider the U(N) gauge group. Note that the U(l) part is decoupled in the 
continuum limit and hence the physics is the same as SU(N) theory. The action is given 
by 



S 



ft K 



i%n—x •En-Eft Vn—yVn—y ynVfi) 2|Xn2/n+:i Vn" 

+ V^2 (oifjaTfjAfj Oifl—x^n x n—x) V^2 (/3n?/nAn fin—y^nyn—y) 
V^2 (On2/n+£Cn ^-n+y^nUrt) V^2 (/3n^n+y^n fin+x£,n%n) 



■n+y\ 



2 2 



2a 2 



+ aV I ynyn 



1 

2^ 



+iVa 2 W a 



V 



N 



1 

2d 2 



i 2 2 

+ a v 



2a 2 



N 



(6) 



Here x, y are N x N complex matrices and a, /3, A, £ are N x N complex Grassmanian 
matrices. They are related to the fields in the continuum by 



1 Xi + iAi 

+ 



4> 



ay/2 
A 



y/2 ' ^ a^/2 
$ = i(a, /3). 



1 X 2 + iA 2 



V2 



(7) 
(8) 



The fermion is in the Weyl representation and is complex. Hence we study the determinant 
of the Dirac operator rather than the Pfaffian. The Dirac operator D is obtained by writing 
the fermion part as ^jiax^ijax^y^H^yi where suffixes i,j, k, I refer to color, a, j3 to spinor 
and x, y to coordinate. 

It is convenient to introduce (semi-)compact decomposition of the bosonic fields [39J 



1 



1 



~ v^T 1 " 1 ' V V2a 
where Ui are unitary, Hi are Hermitian and positive definite, and 

Ui = exp(iaAi), Hi = exp(oXj). 

From x and y, Hi can be obtained as 

Hi = \ / 2a 2 x^x, H 2 = \72aVy. 



(9) 



(10) 



(11) 
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A mass parameter [i gives mass to U(N) scalar fields. The last two terms are not 
present in the original proposal; it gives mass only to U (1) part of the scalar. As observed 
in [191 125] . at large volume and/or with periodic boundary condition for the fermion, flat 
direction in SU(N) sector is dynamically lifted and nonabelian phase (i.e. bound state 
of scalar eigenvalues) becomes meta-stable. (It becomes stabler as N increases.) On the 
other hand U(l) flat direction is never lifted, and in the CKKU model, it can destroy 
the lattice structure. But to stabilize this U(l) flat direction we do not have to turn on 
the U(N) mass (i; the U(l) mass v is fine enough. Given that the U(l) sector is free and 
decouples from the dynamics, nonzero value of v does not affect the supersymmetry i n the 
SU(N) sector in the continuum limit. We explicitly confirm this statement numerically 
(See Figs. El El E|)- Note that this decoupling contains a delicate issue. At lattice 
level, the U(l) and SU(N) secters are not completely decoupled. In order to stabilize 
the lattice structure, the heavy U(l) mass is suitable. However, if the U(l) mass is 
too large, the SUSY breaking effect in the U(l) sector becomes large and it might be 
mediated through the lattice artifact to the SU(N) sector. Our numerical results support 
the decoupling in a wide range of the U(l) mass. It would be nice if the decoupling could 
be explained analytically. We will show that with non-zero value of v the obtained values 
of the observables are the same as those from SU(N) Sugino model (see sec. HJ, which 
justifies the treatment of the U(l) mass term. At finite volume and finite ./V we need 
nonzero /i to completely remove the instability, but at large-A r we can take \x = and 
v ^ so that the supersymmetry in the SU(N) sector is fully restored already at finite 
volume. It enables one to study interesting finite- volume physics like black hole/black 
string phase transition |401 HT] . 

2.2 Structure of light modes 

At the classical vacuum of the U(N) theory, the Dirac operator has 2N 2 fermion zero- 
modes, which correspond to zero-momentum. Apart from the classical vacuum, 2(iV 2 — 1) 
of them are lifted and there remain two zero-modes which correspond to the U(l) part. 
At a discretized level, because of the special property of the CKKU model, only one of 
them is exactly zero [37] • The other approaches to zero in the continuum limit. Let us 
call it the pseudo zero-mode. 

In our simulation, the phase quenched ensemble with det V MMt is used, where M = 
iD. To avoid the exact zero- mode, we add a regulator term to MM^ in the simulation, 

MM f -> MM^ + el. (12) 

In practice e is fixed to be a small enough value (e ~ 1CP 6 ) compared to all nonzero 
eigenvalues of MM^ . 

When we calculate the determinant, we remove the exact zero- mode and the pseudo 
zero- mode by hand. (In practice we remove one exact zero eigenvalue and smallest nonzero 
eigenvalue of the Dirac operator.) As we will see, removal of the pseudo zero-mode is crucial 
to establishing the positivity of the determinant. The reason is simple - in the continuum, 
the Dirac operator has positive determinant because its eigenvalues form pairs (A, A*) after 
removing zero-modes. The reason why we can remove this pseudo zero mode is clear; it 

8 We would like to thank O. Aharony for suggesting the use of 17(1) mass term. 
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will decouple from the dynamics in any case. The corresponding zero-modes which should 
be removed in the lattice simulation are the exact zero-mode and the pseudo zero-mode. 
If the lattice artifact to zero-modes had a pair structure which keeps the positivity of 
the determinant we would not need to remove these two modes. However, phase of the 
pseudo zero-mode fluctuates violently and it dominates the fluctuation of the phase of the 
determinant because pseudo zero-mode does not appear in a pair; without removing it 
the phase of the determinant becomes completely random, just because the phase of the 
pseudo zero-mode is random. 

We numerically calculated the eigenvalues and the determinant of iD. We observed 
that at very small lattice spacing, the other 2(iV 2 — 1) light modes have a pair structure 
(A, —A*), which is exactly expected from the continuum argument. 

2.3 Absence of the sign problem 

As we have explained, it is important to control the flat direction in order to study contin- 
uum physics. In Fig. [I] we show histories of the extent of the scalar fields in the lattice unit 
(left) and in the physical unit (right). When the extent is close to the cutoff scale, 1 in the 
lattice unit, the simulation is not reliable. As we can see from the left panel of Fig. [TJ at 
the physical volume L = 1.0, the U(N) mass fi = 1.0 and the U(l) mass v = 5.0, scalars 
take sufficiently small values, which shows that the flat direction is well under control to 
guarantee the correct continuum limit. Thus there is no instability caused by a lattice 
artifact. Moreover the scalar fields stay finite in physical unit (the right panel). That is, 
this lattice model describes a system without the flat direction in the continuum limit. 
The obtained continuum system has no instability along the flat direction. It is expected 
of course, because we added a mass term. Note that if the scalar converged in the lattice 
unit but diverged in physical unit we would obtain a continuum limit with flat direction, 
which is just a free theory in the Abelian phase. As the mass is decreased, the flat direction 
emerges gradually. In Fig. [2]we show the histories at L = 0.5 with a few values of [i. With 
8 x 8 (a = 0.0625) lattice the fluctuation is not violent even at [i = 0.2 (the left panel). As 
one can easily imagine, with smaller lattice the instability - caused by the lattice artifact 
- appears more easily; see the right panel. This plot uses the same physical parameters 
as the left panel, but a smaller lattice. Note that for the smaller lattice (i.e. larger lattice 
spacing) more spikes appear at fi = 0.2, which is a signal of the instability. The reason is 
obvious - with smaller fi the extent of the scalar in physical unit is larger, and to make it 
to be small in lattice unit lattice spacing must be smaller. By extrapolating to /i = by 
using data from stable region, we obtain a finite extent of the scalar (Fig. [3j here we have 
assumed a simple linear extrapolation, based on the obtained plot). Therefore we can 
expect the phase we are looking at is smoothly connected to the meta-stable non-abelian 
phase [HI [T9l [25] . 

In the following, for all data we show, we have confirmed the simulation does not run 
away to the flat direction. 

Now we are ready to study the sign problem. As discussed in detail in Appendix(XJ the 
determinant of iD has a sign (— l)^ -1 , at least for the constant configurations. Therefore, 
throughout this section, we multiply (— l)^ -1 to the phase factor of the determinant so 
that it localizes around +1. 

In short, what happens both in the CKKU and Sugino is 
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Figure 1: [CKKU] History of the average scalar extent, in the lattice unit (left) and in 
the physical unit (right). At each sweep the average over x and i of ^/Tr(aXi(x)) 2 /N (a 
contraction w.r.t. i is not taken), where aXi = log Hi, is plotted. Clear convergence as 
~ a in the left panel, which corresponds to fixed physical extent of the scalar in the right 
panel, can be seen. The physical volume is L = 1.0, the U(N) mass is /i = 1.0 and the 
U(l) mass is v = 5.0. The gauge group is U{2). 




Figure 2: [CKKU] History of the average scalar extent in the U(2) theory at L = 0.5, 
v = 5.0 and fx = 1.0,0.5 and 0.2. At each sweep the average of \jTr{aXi(x)) 2 /N is 
plotted, (left) 8x8 lattice, even at [i = 0.2 the flat direction is under control, (right) 
6x6 lattice, the flat direction appears around fi = 0.2. Note the difference of the scales 
between the panels. 
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Figure 3: [CKKU] Mass dependence of ^jjTrXf(x)J x (lattice spacing) in the U{2) 
CKKU model. As in Figs. [T] and [21 a contraction w.r.t. i is not taken and an average over 
i and x is taken. The lattice size is 8 x 8, the physical volume is 0.50 x 0.50 and the U{1) 
mass is v = 5.0. 





Figure 4: [CKKU] Distribution of the argument of the determinant in U(2) theory. The 
U(N) mass is \x = 1.0 and the U(l) mass is v = 5.0. The physical volume is L = 0.5 (left) 
and L = 0.75 (right). 



• For fixed lattice spacing, the phase fluctuation becomes smaller at smaller volume 
and/or smaller N. 

• For fixed N and fixed volume, the phase disappears in the continuum limit (small 
lattice spacing). 

Let us start with the U(2) theory. In Fig. H]we have shown how the distributions of the 
argument of the determinant peaks to 0. In each panel, the physical volume is fixed and 
the number of sites is changed. There is a clear tendency that the peak becomes sharper 
as one goes closer to the continuum. 

In order to justify our treatment of the U(l) mass, we have checked the U{1) mass 
dependence of the Wilson loop (Fig. [5|), the norm of the SU(N) part of scalar (Fig. [6J, 
the distribution of the argument of the determinant (Fig. [7J and its cosine (Fig. [S|) . One 
can see that there is almost no U(l) mass dependence. 

In Fig. Owe have plotted the phase distribution at various values of the U(N) mass /x, 
while other parameters are fixed. It turns out that the /i-dependence is small. 
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Figure 5: [CKKU] The 17(1) mass dependence of the Wilson loop in the 17(2) CKKU 
model. The U(N) mass is fixed to \i = 1.0 and the physical volume is fixed to 1.0 x 1.0. 
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Figure 6: [CKKU] The J7(l) mass dependence of (^y jjTrXf(x)) x (lattice spacing) in 

the £7(2) CKKU model. As in Figs. [T]and others, a contraction w.r.t. i is not taken and 
an average over i and x is taken. The U(N) mass is fixed to fi = 1.0 and the physical 
volume is fixed to 1.0 x 1.0. 
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Figure 7: [CKKU] The distribution of the argument of the determinant in the U(2) CKKU 
model, with various values of the U(l) mass v. The left panel is for 4x4 lattice and the 
right is for 8x8 lattice. The U(N) mass is fixed to /i = 1.0 and the physical volume is 
fixed to 1.0 x 1.0. A factor (— l)^ 1 is multiplied. 
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Figure 8: [CKKU] The expectation value of real part of the determinant phase factor vs. 
the U(l) mass u, in the U(2) CKKU model. The U(N) mass is fixed to \i = 1.0 and the 
physical volume is fixed to 1.0 x 1.0. A factor (— l)^ -1 is multiplied. 
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Figure 9: [CKKU] The scalar mass dependence of the argument of the determinant in the 
U(2) theory. The lattice size is 6 x 6, the physical volume is fixed to 0.5 x 0.5. The U(N) 
mass is varied while the U(l) mass is fixed to be v = 5.0. 
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Figure 10: [CKKU] The argument of the determinant in the U(3) theory. The physical 
volume is the L = 1.00, the U(N) mass is \i = 1.0 and the U(l) mass is v = 5.0. 

In Fig. [H)]we plot the phase distribution in the U (3) theory. The distribution is broader 
compared to the U(2) case, but peaks around in the continuum limit. 

In Fig. [11] we plot the real part of the phase. It is clearly seen that it approaches to 
1 as lattice spacing becomes small for each physical volume. The scalar mass dependence 
of the real part is plotted in Fig. [TJl which shows almost no dependence. The detailed 
values are listed in Tables [H [2] and O 

3 Absence of the sign problem in Sugino model 

In this section, we study the Sugino models and observe the argument of the Pfaffian of 
the Dirac operator. As before, we use the scalar mass term to regularize the flat direction 
of the potential. As shown in [25], we have checked that the scalar eigenvalues remain 
close enough to the origin. Therefore we are observing the non-abelian phase, and at the 

9 See Appendix [Bl for the details of the model. 
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Figure 11: [CKKU] The expectation value of the real part of the determinant phase factor 
vs. the lattice spacing. A factor (—) N ~ 1 is multiplied. \x = 1.0, v = 5.0. 
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Figure 12: [CKKU and Sugino] The scalar mass \i dependences of the real part of the phase 
factor in the SU{2) theory (Sugino model, 4x4 lattice) and the U(2) theory (CKKU model, 
6x6 lattice and v = 5.0). The physical volume is 0.5 x 0.5 for both cases. 
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6 
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1.00 


8 
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3 


1.00 


4 


0.748(14) 




1.00 


5 


0.771(20) 




1.00 


6 


0.813(18) 



Table 1: [CKKU] The real part of the phase factor in the U(N) theory. The U(N) mass 
and the U(l) mass are fixed to [i = 1.0 and v = 0.5, respectively. 
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N 


I 1 


Real part of the phase factor 


2 


0.3 


0.943(7) 




0.4 


0.950(6) 




0.5 


0.945(7) 




0.6 


0.958(6) 




0.7 


0.946(8) 




0.8 


0.953(7) 




0.9 


0.965(6) 




1.0 


0.948(8) 



Table 2: [CKKU] The real part of the phase factor in the ^7(2) theory. The physical 
volume is fixed to be 0.5 x 0.5 and the lattice size is 6 x 6. The U(l) mass is fixed to 



N 


V 


Real part of the phase factor 






N x = N y = 4 


N x = N y = 8 


2 


0.0 


0.709(16) 


0.897(10) 




0.2 


0.694(16) 


0.893(10) 




0.4 


0.692(16) 


0.906(09) 




0.6 


0.710(16) 


0.911(08) 




0.8 


0.726(15) 


0.908(09) 




1.0 


0.740(15) 


0.911(08) 




5.0 


0.807(13) 


0.910(09) 




10.0 


0.803(13) 


0.913(09) 




50.0 


0.823(13) 


0.902(10) 



Table 3: [CKKU] The real part of the phase factor in the U(2) theory, with various values 
of U(l) mass v. The physical volume is fixed to be 0.5 x 0.5 and the lattice size is 4 x 4 
and 8x8. The U(N) mass is fixed to fi = 1.0. 
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2-10123 




Figure 13: [Sugino] The argument of the Pfaffian for the SU(2) case. The physical volume 
is fixed to 0.5 x 0.5. The left panel is with various lattice spacings and the right is with 
various scalar masses \i. 



same time, can avoid an unphysical vacuum with large scalar eigenvalues of the cutoff 
scale. The configurations for N = 2 is taken in this work, while those for N = 3, 4, 5 are 
taken from the previous work |25| . 

Let us start with SU(2). In Fig. [131 we plot the distribution of the argument of the 
Pfaffian for fixed physical volume on the left panel. The peak around becomes sharper as 
we go close to the continuum. On the right panel, the scalar mass dependence is plotted. 
Heavier mass gives slightly shaper peak around 0, but the mass-dependence is small. The 
average of the real part of the Pfaffian phase factor is plotted in Fig. [T4l which shows clear 
convergence to 1 as the lattice spacing becomes small. (Note that SU (2) case in the plot 
corresponds to the continuum limit, which uses the fixed volume.) See also Fig. [12]for the 
scalar mass dependence. 

In Fig. [J5J we plot the lattice spacing dependence for the SU (3) theory on the left 
panel. The peak becomes sharper as we go closer to the continuum. On the right panel, 
we show that the phase distribution with fixed lattice size 4x4. The peak becomes sharper 
at smaller volume (or equivalently at smaller lattice spacing). 

The dependence on N is plotted in Fig. [161 As becomes large, the distribution 
spreads. This can be also seen in Figll4| where the lattice spacing dependence of the real 
part of the phase factor is plotted for = 2, 3 and 5. The smaller the lattice spacing is, 
the closer the real part to 1. And larger N shows slower approach to 1. 

The results are listed in terms of the real part of the phase factor in Tables 0] and [5j 

4 Comparison of the CKKU model and the Sugino model 

In order to confirm that our simulation captures the continuum physics, we compare the 
CKKU model and the Sugino model. In the latter, detailed studies have been performed; 
it correctly reproduces analytic results in continuum [25] and also restoration of the full 
supersymmetry has been confirmed [42J. Here we compare the simulation result of the 
CKKU model with the one of the Sugino model with periodic boundary condition for 
fermions. 
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Figure 14: [Sugino] The real part of the phase factor in the 577(2), SU(3) and SU(5) 
theories, with the physical volume 0.236 x 0.236-1.414 x 1.414. 
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Figure 15: [Sugino] The argument of the Pfafhan in the SU(3) theory. The scalar mass 
is fx = 0.20. The left panel is for a fixed volume 0.707 x 0.707 and thus different lattice 
spacings. The right panel is for a fixed 4x4 lattice with various physical volumes (thus 
various lattice spacings). 




Figure 16: [Sugino] The argument of the Pfaffian in the SU(3), SU(A) and SU(5) theories. 
The lattice size is 4 x 4 and the physical volume is 0.707 x 0.707. 
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N 


N X =Ny 


scalar mass [i 


Real part of the phase factor 


2 


2 


0.2 


0.907(4) 




3 


0.2 


0.933(3) 




4 


0.2 


0.953(2) 




4 


0.3 


0.958(2) 




4 


0.4 


0.958(2) 




4 


0.5 


0.968(2) 



Table 4: [Sugino] The real part of the phase factor in the SU(2) theory. The physical 
volume is fixed to be 0.5 x 0.5. 



N 


L 


N X =Ny 


Real part of the phase factor 


3 


0.236 


4 


0.950(3) 




0.354 


4 


0.922(4) 




0.707 


4 


0.857(5) 




0.707 


5 


0.891(14) 




0.707 


6 


0.935(7) 




1.414 


4 


0.742(7) 




1.414 


6 


0.822(14) 


4 


0.707 


4 


0.788(21) 


5 


0.354 


4 


0.82(3) 




0.707 


4 


0.72(4) 




1.414 


4 


0.42(5) 



Table 5: [Sugino] The real part of the phase factor in the SU(3), SU{4) and SU(5) 
theories. The scalar mass is \x = 0.2. 
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Figure 17: [CKKU and Sugino] The expectation value of the Wilson loop (\W\) at [i = 1.0 
and v = 5.0. The extrapolation to the continuum limit has been performed. The gauge 
grope is U{2) and SU(2), respectively. 
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Figure 18: [CKKU and Sugino] The extent of the scalar ^jTr(Xi(x)) 2 /N (contraction w.r.t 
i is not taken). The masses are \i = 1.0 (Sugino and CKKU) and v = 5.0 (CKKU). The 
extrapolation to the continuum limit has been performed. For CKKU, only the SU(N) 
part is plotted. The gauge grope is U(2) and SU(2), respectively. 



Because the simulation of Sugino model was performed with SU (N) gauge group while 
for the CKKU model gauge group was chosen to be U (N) , we compare the absolute value 
of the Wilson loop W = jjTr e l § dxAx : from which the U(l) part decouples. As can be seen 
from Fig. [Tf\ two models give the same result in the continuum limit. We also compare 
the size of the SU (N) part of the scalar fields (Fig. [TBj) . Two models agree reasonably 
well with each other. Therefore we conclude that the both models converge to the same 
continuum limit as expected. 

Before concluding this section, let us comment on the flat direction. In [25] it has been 
shown that the flat direction is lifted dynamically at large- N, both in the continuum theory 
and in the Sugino model, and the Sugino model converges to the correct supersymmetric 
continuum limit. In this section we have seen the CKKU model converges to the same 
limit as well, and hence we can expect the same uplift of the flat direction in the CKKU 
model. 
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5 Why was a "sign problem" observed in previous works? 



In [37] it has been pointed out that the sign of the fermion determinant of N = (2, 2) 
CKKU model fluctuates violently if one chooses randomly generated lattice field configu- 
rations. However, as correctly argued in [57], it does not mean a problem in the continuum 
limit - randomly generated configurations are usually measure zero in the path integral, 
and hence it is necessary to study the distribution of the phase in the phase quenched 
simulation. If the distribution peaks around one (i.e. determinant is real positive) in the 
continuum limit, the sign problem does not exist. More crucial thing is the treatment of 
the pseudo zero- mode. As already mentioned in [37]. removing only the exact zero- mode 
leads to the fluctuation of the phase factor. The pseudo zero-mode, which should give 
zero eigenvalue and should be decoupled in the continuum limit, gives non-zero eigenvalue 
due to the lattice artifact Since what we want to extract from the lattice simulation is 
the continuum limit, we must remove a contribution from the pseudo zero-mode when we 
calculate the determinant. Then we obtain the correct positive determinant. 

Next let us consider a result from an importance sampling in N = (2, 2) theory [381112]. 
According to the plot in the paper, the distribution of the scalar eigenvalues is large; the 
tail of the distribution reaches to (lattice spacing)" 1 . However, in order for the lattice 
considered there to converge to the continuum limit at tree level, the scalar eigenvalues 
must be of order (lattice spacing) . Therefore it is plausible that the simulation does 
not capture the continuum physics 0. In a simulation of the maximally supersymmetric 
theory reported in [38] the distribution of the scalar eigenvalues is narrower compared to 
N = (2, 2) theory, but it is still wide (larger than the size of fluctuation in 4 x 4 lattice in 
Fig- [I]). Therefore it is difficult to obtain robust statements for the continuum limit unless 
studying smaller lattice spacing using bigger lattice. Another subtlety is that in [38] the 
system is projected to SU(N) from U(N). Although it is reported that there is no effect 
to the supersymmetry in this projection, it might have affected to the sign of determinant. 

One interesting observation in [38] for the maximally supersymmetric theory is that 
the phase is close to one when the gauge group is SU(2). The same behavior is observed 
also in one-dimensional theory |23] . In this case there is no apparent kinematic reason like 
the pair structure of eigenvalues, and it is not clear if the absence of sign persists in the 
large volume and/or in the continuum limits. However, if it survives to some extent, it 
may allow detailed study of the SU(2) theory by phase-quench or reweighting. 

In = (4, 4) theory, it has been reported that a lattice model a la CKKU suffers 
from the sign problem [44] . This 8 supercharge system is known to have the sign problem 
in general so the result itself is reasonable. It is interesting to see whether the result 
changes when pseudo zero-mode is removed, although we do not expect the pair structure 
of eigenvalues in this theory. In addition, it is not clear whether the flat direction was 
under control in [44] . which is an important point to be studied. 



As discussed in |37| . one can regard this artifact is caused by the orbifold projection, which does not 
commute with manipulation needed to prove the positivity of the determinant. 

11 It has been remarked that large phase fluctuation arise when scalars take large expectation values 
[43], and that this "fluctuation of sign" suggests the SUSY breaking because it can make the Witten index 
vanish. However, if such configuration corresponded to continuum theory, it must be an abelian phase, 
which does not have any dynamics, and hence the supersymmetry cannot be broken. 
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6 Conclusion 



In this paper we established the absence of the sign problem in 2d N = (2, 2) super Yang- 
Mills theory on lattice. We studied two different lattice models: The Sugino model and 
the Cohen-Kaplan-Katz-Unsal (CKKU) model. We have clarified the structure of the light 
modes in the CKKU model and pointed out the importance of the removal of the pseudo 
zero mode. We also confirmed that the both lattice models provide the same continuum 
physics as expected. 

As we pointed out in § [5j and also discussed in [25J, in order to obtain correct continuum 
limit it is crucial to control the scalar flat direction. From this point of view, it is possible 
that some of the past simulations for two-dimensional super Yang-Mills theories failed to 
capture the continuum physics. It is urgent to check whether the scalar flat direction was 
under control. 
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A An overall sign of the determinant of Dirac operator in 
the CKKU model 

Let us introduce the following notation, 

Vf{n) = V2atv(T a x H T b ), (13) 

Vf{n) = V2atr(T a x H T b ), (14) 

V 2 ab (n) = V2atv(T a y H T b ), (15) 

V 2 ab (n) = V2atv{T a y n T b ), (16) 

where T a is a Hermitian gauge generator normalized as tr(T a T b ) = 5 ab . and are 
related by 

% = (yp T = 4 (17) 

The fermionic part of the action is expressed aJ^l 

S F = i(a a H , /3S) (?m\<* ?m\<*) (f) = i(a a H , ($) (18) 



rjab 77 I \ c I v n.' ■ n) nm \ c 

nm nm m/ W 



A factor N/a in front of D, which is irrelevant to the argument here, is omitted. 
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with 

Dfrnkx = ~i (S^rnMn) ~ (n))^ , (19) 

Dmkx = -i (S H , A V 2 (H) - S^V 2 *(n)y b , (20) 

DfrnU = ~i (-5ft,rnV 2 (fi + 1) + 6~^V 2 *(n - 2))°" , (21) 

Dm\M = ~i (Sn,Mn + 2) - ^^VT (n - 1))^ , (22) 
where each of the fermion components is defined as 

a n = T a a% etc. (23) 
Assuming the bosonic fields are constants we obtain the momentum representation as 

Dtp) = -i ( Vl ~ ^ V * ~ V * + e: iaP2 V 2 *\ (24 ) 
We further decompose it into 



D (p) = \ 2 ?a P2 ) D'(p) [ «^ MP2 ] , (25) 

\ e — —) \0 e 2 + 2 J 

where D' has the same determinant as D. The explicit form of D' is 

/ iap-\ iap-\ iapn iapn \ 

D'(v) - -i V 1 ~ ~ V ^ ~ V * ~ + M (26) 

\Jr ) I iap22 ^ a P2 iap-^ iap^ I • J 



V 2 e — — - V*e— Vi e~ - Vfe 



2 



Note that because of a factor a/2 in the exponents, D'(p) is not periodic w.r.t. 2ir/a so 
that we have to be careful about the treatment of the the boundary of the Brillouin zone. 
We use a region —ir/a < p-i < it /a, where — -n/a is not included. It is easy to see that it 
satisfies 

a 2 D'(p)a 2 = D'*(-p). (27) 

Therefore, for p which has — p in our Brillouin zone, we have the following quartet of the 
eigenvalues P*l 

\(p),\(-p),\*(p),\*(-p) (28) 

and the contribution to the determinant is always positive. For p = 0, we have a pair 
A(0), A*(0) which has a positive contribution as well. Note that for p = 0, D(p) = D'(p) 
thus the eigenvalues of D(p = 0) make a pair (A, A*) as well. 

If the lattice is (odd) x (odd) , we do not have pi = it /a modes so the determinant of 
D' and thus that of D is positive. 

If the lattice is (even) x (even), since we have pi = it /a modes and p% = —ir/a is not 
in the Brillouin zone, the sign becomes non-trivial. For p\ = n/a, we have 

a 2 D'(ir/a,p 2 )a 2 = D'(ir/a, -p 2 )* (29) 



Do not confuse eigenvalue A with one of the fermion or 't Hooft coupling. 
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for p2 7^ 7r/2. Therefore we have a quartet (A(p2), A(— P2); A(p2)*> A*(— P2)) for P2 / and 
a doublet (A, A*) for p 2 = 0, both of them give positive contributions to the determinant. 
The situation is the same when P2 = ir/a. The remaining combination is p = (ir/a,Tr/a). 
This time we have 

D (p! = tt/o, p 2 = vr/a) = " + ^ _^ _ ^ ) (30) 

which satisfies U2D' a% = —D'*. The eigenvalues make a pair (A, —A*). The sign contribu- 
tion to the determinant from this sector is (— 1)^ = (— 1)^. 

For iD, an extra factor ^ 2Ar2x ( num - of lattice sites) a pp earS) anc i hence the phase is (— 1)^ 
for both odd x odd and even x even lattices. Since we remove the two lightest modes - 
one exact zero mode and one pseudo zero mode — the determinant of iD picks up an 
extra factor i~ 2 and thus the determinant of iD has a sign (— l)^ -1 . 

B The action of the Sugino model 

Sugino's lattice action [8] is given by 



Siattice = a x a y I Cbi{x) + £Fi(£) f + (auxiliary field), 

x \ i=l i=l ) 



(31) 



where 



C B i{x) 


N 


S Gl rj* Qj^j 




N 


C-B2{x) 


8a 2 a 2 




N 


2a%a y 



Tr[0(x),0(x)] 2 , (32) 

Tr$ TL (x)\ (33) 

Tr{ (4>(x) - U x (x)(f>(x + a x x)U x {x)~ 1 ) 

x (4>(x) - U x {x)4>{x + a^U^x)" 1 ) } 
+ ^3 Tr i(^ " Uyix^x + ayy)Uy(x)- 1 ) 

ZCL x CLy 

x - U y (x)4>{x + a^y)^^)- 1 ) } (34) 



14 Here we follow the notation in |32| with a slightly different normalization. Although group theoretical 
normalizations are different in [32] and [7], one can absorb them by rescaling the 't Hooft couplings as 
2Aqkku = Asugino- We set A = Ackku = 1. 
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and 

N 

C F i(x) = -—^Tr(rj(x)[<j>(x),ri(x))), (35) 
N 

Cf2{x) = - 9 2 2 Tr(x{x)[(/)(x),x(x)]) , (36) 
£-F3(x) = - r-g- Tr (^o(x)V'o(^) (<ft(s) + U x {x)4>{x + a^L^x)" 1 ) } , (37) 

C F4 (x) = --^Tr{Mx)Mx)(Hx) + U y (x)4>(x + a y y)U y {x)- 1 )} , (38) 

N / \ 
£fs(x) = i^2 Tr [X(x) • Q$(x)J , (39) 

£f6(£) = y Tr{^o - U x (x)r](x + a^x ^(x) -1 ) } 

-t— gTr{^i (77(a) - U y (x)ri(x + ayy)Uy(x)- 1 ) }, (40) 

where U(x,fi) are gauge link variables, </>(x) is a complex scalar, 77(a), x(^) an d ipu,{x) are 
fermion field, and are lattice spacings El, e is a real parameter which must be chosen 
appropriately for each N, 

* (5 ° = ~l-\f-PWW* ' ^) = ^)-^( T ^))- 1 ' W 

where -P(cc) = U x (x)U y (x + x)ul{x + y)Uy(x) is the plaquette variable, and Q generates 
one of the four super transformations, 



QU^x) 


= iip^fflU^x), 


(42) 


Qi>n(x) 


= iippffltp^x) - i(4>(x) - U^(x)4>(x + a M /l)C/ M (x) _1 ), 


(43) 


Qm 


= 0, 


(44) 


Qx(x) 


= H(x), 


(45) 


QH{x) 


= [</>(x),x(?)]> 


(46) 


Q4>(x) 


= T)(x), 


(47) 


Qrj(x) 


= [</>(£) J(x)]. 


(48) 



Sugino's action Sluice is invariant under the supersymmetry generated by Q, because Q 
is nilpotent up to commutators and S can be written in a Q-exact form. 

In [8], using super-renormalizability and symmetry argument, it was shown that other 
three supersymmetries, which is broken by a lattice artifact at the discretized level, is 
restored in the continuum limit. Furthermore, in |42j . this restoration has been confirmed 
explicitly by the Monte-Carlo simulation. Absence of operator mixing/renormalization is 
has been shown perturbatively in 



15 In the actual simulation we have used the isotropic lattice, a x 
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